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ABSTRACT 

V-&I  this  paper**  wo  contider&the  lot  size  model  for  die  prodnedoo  and  storage 
of  a  tingle  commodity  with  timitatiooi  on  production  capacity  and  the 
peatibiHty  of  not  meeting  demand,  Le.,  stockouts,  at  a  penalty.  The  stockout 
option  means  that  horizons  can  exist  and  permits  the  use  of  horizons  to 
develop  ft  forward  algorithm  for  solving  the  problem.  The  forward  algorithm 
ir  shown  in  the  wont  caae  to  be  asymptotically  linear  in  computational 
wfihmgmta,  in  contrast  to  the  case  for  the  clestical  lot  tize  model  which  has 
eatponential  computing  requirements.  Two  versions  of  the  model  are 
considered:  the  first  in  which  die  upper  bound  on  production  is  the  —me  for 
ovary  time  period;  and  the  second,  in  which  the  upper  bound  on  production 
is  permitted  to  vary  each  time  period.  In  the  first  case  the  wont  case 
computational  difficulty  inert—  in  a  cubic  fashion  initially,  and  then 
beooeam  linear.  In  the  aecoad  earn  the  initial  increase  is  eatponential  before 
becoming  timer.  Batides  the  forward  algorithm,  a  number  of  mnnmatj 
conditions  an  derived  which  reduce  the  computational  burden  of  solving  the 
iniagsr  programming  problem  posed  by  the  model. 
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L  INTRODUCTION 

The  elamie  lot  size  model  [93  involves  the  production  of  a  single  product, 
merit  it  hi  a  warehoura  of  unlimited  capacity  until  needed,  and  require*  the  complete 
futfQhnent  of  all  (detenniniatic)  demand*.  Various  modifications  have  been  made  to 
that  model  which  include  the  introduction  of  upper  bounds  cm  production  and  on 
inventory  apace,  and  the  w>tngynt  of  orders.  Each  modification  permits  the 
derivation  of  special  results  relevant  to  that  case.  In  die  present  paper  we  consider  a 
new  variation  of  the  model  that  of  permitting  stockouts. 

Owe  Has  of  rswerch  in  the  area  of  production  planning  deals  with  capacitated 
lot  sine  rtefsBi  of  a  single  item.  Once  the  single  item  case  is  understood, 
kaowiadgs  from  it  is  used  to  solve  the  multiple  item  situation.  Research  in  this  area 
was  ptonswad  by  the  uncapadtarad  lot  sine  modal  of  Wagner  and  Whitin  [93.  Others 
hMhida  Zabot  [103.  Z**wtH  [113  [123.  and  Lundin  sad  Morion  [73.  Each  of  theae 
rafrrtoas  iavotvm  a  forward  algorithm  approach  which  aolvw  succmsively  loepr  finite 
horizon  probisma  until  a  decision  horizon  is  snoounrarad.  By  definition,  if  optimal 
p» adoption  derisions  daring  dm  interval  [U*3  are  completely  independent*  of  the 
daasand  drae  bayoad  t**  h  t*.  than  t*  is  a  daririea  horizon  and  t**  is  a  forecast 
horizon.  Ptem  a  practical  standpoint,  wa  are  really  only  inter  rated  in  determining  the 
fint  daeMoo  horizon  abaca  the  production  problem  is  typically  resolved  periodically  to 
Incorporate  improved  forecast  data.  Single  item  capacitated  lot  size  model  research  has 
bean  Untiled  to  static  horizon  length  probisms  and  raramptioas  as  to  the  demand 
pattern,  capacity  pattern,  rad  cost  functions  (era  Florian  and  Klein  [23.  Jagnmefhen 
and  Rao  [33.  Love  [f3.  Swovehmd  [S3.  Louveenx  [53.  Korpooker  [43.  and  Baker. 
Dixon,  Magazine,  and  Silver  [13).  The  drawback  of  the  static  horizon  assumption  is 
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that  information  from  beyond  the  ead  of  the  horizon  could  chant  the  entire 
prwdnrtkm  plan;  that  k,  oo  dadrion  horizoes  are  obtainable 

The  purpcac  of  this  paper  is  to  tow  how  decWon  and  forecast  horizons  nay 
be  obtained  for  single  item  capacitated  lot  size  models  In  order  to  obtain  horizons, 
atocfcouts  dost  sales)  are  permitted  but  no  backorders.  This  situation  has  not  bean 
tionsidared  by  previous  authors.  A  concave  production  function  nonsisting  of  e  find 
coat  component  end  e  How  component  is  assumed.  Production  is  bounded  by  capacity 

OOBBifUnll  WAK41  ul*7  Or  IQftj  DOl  06  COQSWU  irOB  pOnOO  ID  pMKXL 

The  formulation  of  the  modal  under  study  is  pnmntsd  in  Section  1  Mar  emery 
conditions  for  a  solution  to  be  optimal  are  stated  and  proved  in  Section  3.  Section  4 
two  types  of  horizon  theorems  are  discussed.  The  first  type  involves  dm  interaction  of 
the  tneWwt  and  hoMfaf  cotta,  the  ti—wrf  axiom  booms  of  «*—*»■*■»«»  on  production 
capacity.  Two  forward  solution  algorithms  that  uriHat  both  dm  optimality  conditions 
and  dm  horizon  theorems  are  prammd  hi  Section  5.  The  first  apptias  to  problems  for 
whmi  un  guumium  pfocimKio  ■  nonii  in m  pwwo  10  pmoo.  11  m  nown 

that  tba  wont  can  effort  involved  in  solving  this  modal  is  initially  cubic  but 
mympeotically  baconrna  linear  in  dm  length  of  dm  problwn  The  second  forward 
upnan  trma  xm  cm  to  whjcii  qw  mssmumi  prooucuo#  cip>ciiy  hm  wiia  imm, 
and  bare  it  is  found  that  the  wont  can  effort  is  initially  exponential  but  asymptotically 
bcooams  linear  in  the  length  of  the  problem.  If  decision  horizons  are  detected 
frequently,  so  that  dm  lengths  of  dm  snbproMems  are  kept  small  even  this  problem  is 
not  too  computationally  demanding.  Computational  results  with  the  algorithms  ere  given 
in  Section  i. 
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1  MODEL  FORMULATION 

The  capecitned  tot  size  model  with  stockouts,  which  is  ta  intapr  liaeer 
program,  is  glean  hi  Figure  2-1  and  is  called  Problem  PT.  la  Problem  PT,  Objective 
faacttoa  1  is  to  be  mjahaiaed  subject  to  Coostratots  (2)  -  (5). 

Mia  Z*M<#4r  ♦  fXr  *  hlr  ♦  *r)  (1) 


IrH  ♦  Xr  -  Ir  ♦  Sf  *  dr  for  r-U..T  (2) 

Xr  -  XMAXr4r  *  0  for  r»L-..T  (3) 

Xr,  Ir.  Sr  noimepttvs  sad  integer  for  r“L~..T  (4) 

binary  for  r*L.„,T  (5) 


Figure  2-1:  Problem  PT:  Capacitated  Lot  Size  Model  with  Stockouts 

The  eariabtoa  of  the  model  are:  Xf.  the  productioa  quaatity  in  period  r\  1^; 
the  Inventory  toed  at  the  ead  of  period  r;  Sr>  the  number  of  stockouts  incurred  in 
period  r;  and  if,  a  binary  variable  which  is  1  when  Xf  >  0  and  0  otherwise;  The 
parameters  of  the  model  are:  the  setup  cost;  p.  the  variable  production  cost  per  unit 
(assumed  constant);  h.  the  holding  cost  per  unit  per  period;  s.  the  stockout  cost  per 
unit;  XMAXr,  the  production  capacity  in  period  r;  and  the  demand  requirements 
ia  period  r. 
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The  objective  it  to  minimize  the  sum  of  total  aetup  coats,  total  variable 

yntetin  casts.  total  cotta,  and  total  stockout  costs  as  stated  in  Objective 

t 

Puncthm  (D.  Constraint  (2)  defines  the  production-invantory-stockout  relationship  with 
the  daooand  raquimasota  It  is  seamed  that  I  ■  I  *  0  without  loot  of  generality  and 
that  Ir  fe  0  for  all  t  so  that  backorders  are  prohibited.  The  upper  boundfcj  erf 
production  is  scoompliahed  by  constraint  (3).  Constraint  (3)  also  forces  the  setup  cost  to 
be  incurred  whan  production  is  positive.  Constraint  (4)  requires  the  production, 
inventory,  and  stockout  quantities  in  each  period  to  be  noonegative  and  integer. 
Finally,  constraint  (5)  impoass  the  binary  {0.11  restriction  on  the  S  variables. 

Probtan  PT  can  be  formulated  as  a  concave  cost  network  flow  problem  as 
follows.  Let  node  M  be  the  master  supply  node  that  contains  the  sum  of  die  supplies 
for  all  T  periods.  We  define  two  transshipment  nodes.  P  and  S.  Node  P  will  transship 
aS  the  units  that  mtfefy  demenda  via  production  and  node  S  will  tranship  all  the  units 
that  satisfy  dsnmnds  via  stockouta  Finally,  nodes  1.  2.  ....  T  ere  die  period  demand 
nodes  such  that  the  demand  at  node  r  is  &r-  The  arc  set  consists  of  the  following 
four  are  subasttc 

L  Directed  arcs  (M.P)  and  (MJS)  which  ship  flows  at  zero  cost  and  have  upper 
bounds  of  «oo: 

2.  Directed  arcs  (P.r)  with  upper  bounds  of  XMAXr  for  r  *  1 . T 

which  each  incur  zero  shipping  cost  for  zero  arc  flow  and  e  cost  of  a  ♦  (p 
•  Xr)  if  die  arc’s  flow.  Xr>  is  greater  than  zero  (then  constitute  the 
production  arcs); 

3.  Directed  arcs  (S,r)  with  upper  bounds  of  *oo  for  r  *  1 . T  which 


at tip  flows  at  a  cost  of  s  per  unit  flow  (these  are  the  stockout  arcs}; 


4  Directed  arcs  (r,r+D  with  upper  bounds  of  »  for  r  »  1.  .  .  .  .  T-l  which 
ship  flows  at  a  cost  ot  h  pa  unit  flow  (the  inventory  arcs). 

A  pictorial  representation  of  the  Problem  PT  flow  network  is  given  in  Figure  2*1 


An  optimal  solution  to  Problem  PT  will  be  a  spanning  tree  of  the  PT  network 
in  Figure  2*2  and  each  unit  travels  from  the  supply  node,  M,  to  its  respective  demand 
node  by  the  lowest  cost  available  path. 


3.  NECESSARY  CONDITIONS  FOR  A  SOLUTION  TO  BE 
OPTIMAL 

Nsc— ry  conditions  on  the  optimal  solution  to  Problem  PT  are  now  stated  and 

proved. 

Definition  1:  Define  die  maximum  holding  pariod  to  be  k  = 
L(s-p)/bJ.  The  following  result  justifies  tills  name. 

2:  In  an  optimal  solution,  the  marinium  number  of  periods  a 
unit  will  be  held  in  inventory  is  the  maximum  holding  period  k. 

Proof:  Assume  there  exists  an  optimal  solution  with  one  or  more 
units  being  held  in  inventory  for  more  than  k  periods,  say  k  +  1  periods. 
Than  there  will  be  a  path  of  k  +  1  basic  inventory  arcs  from  a  period  t  to 
pariod  t  *  k  ♦  L  The  total  bolding  and  variable  production  cost  for  one  unit 
will  be  p  ♦  h*(k*l).  Now  since  k  »  L(s-p)/hJ.  which  implies  t(s-p)/h]  -  l 
<  k  S  (s-p)/h,  or.  solving  for  a, 

p  ♦  (h*k)  S$<p  +  h*(k+l). 

The  rightmost  inequality  implies  that  the  stockout  cost  is  less  than  the  sum  of 
its  total  holding  cost  and  variable  production  cost  Therefore  the  original 
solution  could  not  have  been  optimal.  ■ 

The  lemma  which  follows  describes  how  units  should  be  distributed  from 
inventory  in  an  optimal  solution  to  Problem  PT.  The  FIFO  rule,  or  first-in.  first-out 
rule,  distributes  the  oldest  item  first  while  the  LIFO  rule  (last-in,  first-out  rule) 
distributes  the  newest  item  first  Both  the  FIFO  and  UFO  rules  should  be  regarded  as 
bookkeeping  systems  for  inventory  valuation.  Rarely  are  these  systems  used  as  the 
actual  method  of  distributing  units  from  inventory.  The  main  purpose  for  proving  this 
lemma  is  to  aid  in  the  proof  of  an  important  theorem  which  follows  it 


3:  la  an  optimal  solution,  the  newest  item  in  inventory 


should  be  distributed  first  (LIFO). 

Proof:  We  begin  with  a  plan  that  for  at  least  two  units,  the  first 
unit  produced  is  distributed  before  the  second  (FIFO).  It  is  then  shown  that 
switching  the  distribution  of  these  two  units  from  FIFO  to  LIFO  will  never 
lend  to  a  higher  cost  solution. 


Consider  two  units  that  are  distributed  on  a  FIFO  basis: 

Unit  When  Produced  When  Distributed 


It  is  smuined  that  £  t}  (which  implies  that  tj  S  t*  for  FIFO  to 

bold)  and  that  t’  -  tj  h  k  and  tj  -  t2  *  k  so  that  both  units  ere  produced. 
Case  is  tt  *  tf  or  t*  *  t* 

Case  II:  tt  *  tf  *  t’  «  t’ 

Case  Ilia:  t,  +  k  h  t#' 

Case  Illb:  t,  +  k  <  t#' 

For  Cases  I  and  II,  the  total  costs  under  FIFO  and  LIFO  are  equal 


For  Case  Ilia  and  Illb,  t  <  t2  £  t|  <  t\  The  FIFO  total  variable 
production  and  holding  cost  for  the  two  units  is:  2*p  +  h*C(t*-t  )3. 

If  we  switch  to  LIFO,  the  distribution  pattern  will  be  as  follows: 

Unit  When  Produced  When  Distributed 

If  tf  +  k  h  t’  (Case  IHa),  then  unit  i(  will  still  be  produced  under 
UFO.  The  LIFO  total  variable  production  and  holding  cost  for  the  two  units 


k  2*p  +  b*[(t*-t  Mlj-tj)]  which  is  equal  to  the  FIFO  boldine  cost, 

If  t(  ♦  k  <  t’  (Case  mb),  then  unit  if  will  not  be  produced  under 
LIFO.  The  LIFO  total  variable  production,  holding,  and  stockout  cost  is  p  ♦ 
h*(tj-t )  ♦  s.  The  proof  is  complete  when  it  is  shown  that  p  +  h*(f-t )  +  s 
S  2*p  ♦  h+KtMWr-yi. 

Assume  not,  Le..  p  +  h*(tj-ta)  +  s  >  2*p  +  h*[(t*-t  Mf£-t  )3,  which 
reduces  to: 

h  *  (tVt)  <  s  -  p  (6) 

By  definition,  k  *  L(s-p)/hJ,  which  implies  that 

h*ki  s-p<h*  (k+1).  (7) 

Recall  that  it  was  assumed  that  t(  +  k  <  t*  or 

k  <  t*  -  t .  (8) 

Combining  inequalities  (6)  and  (7)  yields  h  •  (t’-t)  <  h  *  (k+1)  and 
since  it  was  assumed  that  h  >  0, 

r2  -  t,  <  k  +  L  (9) 

Combining  inequalities  8  and  9  leads  to 
k<t’-t)<k  +  l 

which  is  impossible  since  k.  t',  and  t(  are  integers.  ■ 

Theorem  4:  In  an  optimal  solution,  every  k  +  1  consecutive  periods 
forms  a  regeneration  set,  Le.,  contains  one  or  more  regeneration  points  at 
which  the  ending  inventory  is  zero. 


Proof:  Assume  the  contrary,  that  there  exists  an  optimal  solution  with 
k  ♦  1  conm eutfre  periods  that  do  not  form  a  regeneration  set,  say  periods  t’ 
-k.t*-k«-l.....T  +  l.  so  that  there  are  k  +  1  periods  of  positive 
endtag  inventory.  Since  an  optimal  plan  can  be  found  using  LIFO 
distribution  (Lemma  3),  a  unit  produced  in  period  t’  -  k  will  be  held  until  at 
least  period  V  *  L  that  is,  at  least  k  ♦  1  periods,  which  contradicts  Lemma 
1  • 

Theorem  4  states  that  in  an  optimal  solution,  for  a  sequence  of  any  k  +  1 
consecutive  periods,  the  inventory  level  at  the  end  of  at  least  one  of  these  periods  must 
be  aero  (Le.,  a  regeneration  point).  The  determination  of  a  regeneration  point  makes  it 
possible  to  partition  a  large  problem  into  smaller  subproblems,  each  of  which  can  be 
solved  easily. 

lemma  &  In  an  optimal  solution,  1^  •  S  *  0  for  all  L 

Proof:  If  in  an  optimal  solution.  I  ,  •  S  ,  >  0  for  some  t\  then  there 

is  a  path  from  node  M  to  node  t  ♦  1  via  nodes  S  and  t  at  cost  s  +  h.  The 

path  via  node  S  only,  using  arc  (S.t+1),  has  cost  s,  so  that  a  lower  cost 

solution  is  obtainable  by  using  this  alternate  route,  contradicting  the 

assumption  that  the  original  solution  was  optimal.  ■ 

Definition  6:  Define  the  minimum  positive  flow  on  a  production 
arc  to  be  XMIN  *  [^/(s-p)],  where  M  is  defined  to  be  the  smallest  integer 
>*  x.  The  next  theorem  justifies  this  name. 

Theorem  7:  In  an  optimal  solution,  if  X  >  0,  then  XMIN  £  X  £ 

l  t 

XMAX. 

t 

Proof:  Gearly  X  £  XMAX  for  all  L  Now  assume  there  is  an 
optimal  solution  such  that  on  a  production  arc  (P.t*)  there  is  a  positive  flow 


X.  <  W(a-p).  (The  oeiUns  brackets  nay  be  removed  tines  X  most  be 

*  •  I 

integnr.)  Rearranging.  we  have  s  •  X,  <  9  +  (p*X ,).  This  means  it  is 
cheaper  to  stock  out  of  the  X,  units  than  it  is  to  produce  them,  a 
contradiction.  ■ 

Corollary  S:  If  XMAX,  <  XMIN  for  some  f,  then  X(,  *  0. 

A  theorem  similar  to  the  next  one  is  proved  in  Baker.  Dixon.  Magazine,  and 
Silver  [1]. 

Theorem  9:  In  an  optimal  solution.  I  •  (XMAX-X)  *  X  -  0  for 

aD  t 

Proof:  Assume  die  contrary,  diet  there  exists  a  period  t*  in  an 
optimal  solution  with  I.  *  (XMAX-X.)  •  X.  >  0.  This  means  that  arcs 
(f— l*f)  and  (P.f)  are  both  basic  (refer  to  Figure  2-2).  Since  arc  (f-l,t*)  is 
basic,  there  must  be  an  arc  (P.r)  (t"<f)  with  a  positive  flow.  This  positive 
flow  must  be  at  XMAX„,  otherwise  a  cycle  win  have  formed,  violating  the 
characterization  of  optimal  flows  as  ying  trees.  However,  a  lower  cost 
flow  is  obtainable  by  increating  the  flow  on  (P.f)  until  it  becomes  noobatic 
and  decreasing  the  flow  on  (P.f*)  by  the  amount  of  die  increase  on  (P.f). 
Production  costs  remain  the  same  but  total  holding  costs  are  reduced.  Since  a 
lower  cost  flow  is  poatible,  this  contradicts  the  initial  assumption.  ■ 
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Itaaaaoa  ttha  ifW  nMh,  if  k  -  l(a-p»/kj  <  la-pf/h  aai 
tar  mm  p atai  L  I  >  I  tai  f«  wl  r  h  ta  Mortal  H  1  r  1  A  if 
X  >  •  taaa  Xr  *  XMAX,. 

*  >  ft  m«  a  pariai  tr  (r-k  irirl* kfc  ft  <  X_  <  XMAX  . 

r  r 


v 


•  (r\r**a.  <r-*ur*B. . 
r  -  r*.  titao  a  epda 


,  (r-LD  am  afi 

ia 


Now  mhm  that  at  laaai  aoa  of  taa  ternary  area  tamo  porta* 
r*  aod  r  <f  <  r)  ia  aaohNte  vita  am  flaw.  Soar  *r  >  0  aai  ft  <  X  „ 
<  XMAX^  it  aowt  ta  etaapar  la  aock  a*  la  pariai  r  * ■  v  teenaaa 
productfoo  io  pariai  r*  ta  aatafy  iaaMota  io  pariai  r.  Ttat  ta 

p  ♦  *»<r-r*)  >  a 

Now  rioea  k  -  l(*-p»/hJ  <  (t-pi/k  arir-kir.ww  otaata 

p  ♦  taCr-tr)  <  a 


which  ia  a  oaouaiictioo.  Ttaraforo  oar  origtaal  aotarioo  oooM  oat 


Tta  aaxt  thaorwo  ia  a  variation  oa  ooo  doa  to  Bakar,  Dixoo,  Maprtoa, 


Stfror  £13 


lit  to  a 


X>  •  mBUMAX,,  *  not  pra 
mxr- aaoEMAX^.  3^4)  *  xmin. 


1 1 1  x  >  01. 


r  ♦  k  s  t.  if  r  ♦  k  >  T. 


Let  f  ■  m<  t  I  X  >  I  },  U.  r  k  tbs  a  period  of 


ad  nano  dnt  f  ♦  I  i  T.  (The  proof  cea  be  eerily 


arid  idea  an  r  ♦  k  >  TL  Oeariy  XMAX .  iX  H  If 
XMAX',  <  Xj^d,  riM  X,  an  not  a  X^dj  ad  period  r  ♦  k  mwt  bare 


XMAX.  this 


MM  Ml  tf  Xpdj  <  XMAX,.  dai  astd*  X.  equal  to 
riai  Xf*d,  wB  *ve  a  Mtaria  avtag  «  at  k  ♦  1  periods 

K  j 


ria  aria  *t  at  k  ♦  1  periods 

4  If  X.  is  a  to  some  quantity 


te  parted  t*  ♦  k.  which 


a  ■ 


ix  a 


XMAX  a«  X  >  4  riM  X  -  XMAX. 

•  i  ii 


if  for 


parted  t,  d  at 


ctmbx  Mm  mv  no  io  opmi  wroon  ditB|  ft  pftnoa  r 

wMh  dr  k  XMAX,  ad  0  <  X.  <  XMAX..  Tan.  either  area  (t*-l.f)  and 


CP.f)  an  aili  aria  (wtech  contradicts 


9),  or  am  (S.t*)  and  CP.t’) 


10).  Heoce  the  theorem  is  true,  a 


aaa  above,  it  is  possible  to 


conditions  an  optimal 


aaat  satisfy.  A  solution  that  ntisfa  these  conditions  is  not  necessarily 
iA,  than  conditions  are  not  sufficient  for  optimality.  However  use  of  the 
'  conditions  redeem  die  aim  of  the  solution  space  diet  must  a  considered  in 


for  a  optimal  solution. 
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4  THE  DETECTION  OF  DECISION  AND  FORECAST 
HORIZONS 

In  this  section  throe  theorems  are  preaented  that  sometimes  permit  the 
ftejermtearinn  of  dedrion  and  forecast  horizons  for  the  problem  under  consideration. 
Whan  they  exist,  them  horizons  are  used  in  the  forward  algorithms  described  in  Section 
5. 

A  forward  algorithm  begins  by  solving  die  initial  one  period  problem  and 
proceeds  by  solving  imck— ivaly  longer  finite  horizon  problems.  To  be  useful  there 
moat  be  a  mopping  rule  in  order  to  know  when  to  terminate  the  forward  algorithm. 
Whan  decision  and  forecast  horizons  can  be  detected  die  algorithm  can  stop,  since  the 
solution  up  to  die  decision  horizon  is  guaranteed  to  be  pert  of  the  optimal  solution  to 
any  problem  longer  than  the  forecast  horizon.  Specifically,  the  solution  for  the  first 
period  win  be  pert  of  the  optimal  solution  to  any  problem  longer  than  the  decision 
horizon.  Sian  production  planning  models  are  unally  re-solved  every  few  periods  to 
make  uw  of  updated,  more  accurate  forecast  information,  only  the  first  few  decisions 
are  needed. 

Deflnitiea  13s  If  Xj  is  the  optimal  period  j  production  quantity  in 
the  t-period  problem  P\  than  the  sequenn  1X‘,  ....  X1^  is  called  an 
optimal  production  sequence  for  Problem  P*. 

Definition  14:  If  the  sequenn  (X*.  .  .  .  .  Xjl  is  an  optimal' 
production  sequenn  for  Problem  P‘,  then  the  sequenn  IX*.  ....  X‘,1  is  an 
optimal  production  aubcoquonco  for  Problem  P1,  provided  t’  £  t 

Theorem  IS:  Consider  problems  P1* . P1  where  k  is  the 

maximum  holding  period  (k  ♦  1  £  t  £  T).  If  Xp  ■  Xp*‘  *  ...»  X1  for 


Pm ft  Prom  Theorem  4,  it  k  known  teat  for  any  aaqoaaaa  of  k  ♦  1 


oooaeeutm  p^kdt  tba  ending  inventory  of  at  team  ooa  of  theee  periods 
mum  ba  zero  in  an  optimal  rotation.  Whan  rotates  protean*  P*"*,  F***1,  .  .  . 

,  P‘,  tee  inventory  at  tee  end  of  each  proteam  hat  to  be  mo  in  any  optimal 
solution.  Therefore,  the  production  sequence  in  at  team  ooa  of  thaae 
problems  will  be  part  of  die  optimal  solution  to  any  problem  of  length  t*. 
where  t  £  t*  S  T.  If  X*  «  Xj**'  *...«Xjforj»l.....rtt*r 

£  t  -  k),  this  means  that  each  problem  P*"* . P*  has  the  same  optimal 

production  subeequence  for  the  firm  f  periods  and  therefore  this  subeaquaeoe 
mum  be  part  of  die  optimal  solution  to  any  problem  of  length  t\  where  t  1 
s  i.  lnereiore  pcnoG  r  is  t  osctaon  uonzocL  atoot  i  ponoos  01  omudri 
tn  nasded  to  dstsnniM  ^  dacWon  bofino,  pariod  I  is  a 
forecam  horizon,  a 

This  theorem  te  similar  to  the  Lundta-Morton  decision  horizon  theorem  [7]. 
We  know  from  Theorem  4  diet  every  k  ♦  1  consecutive  periods  form  a  repnsrifion 
sat  If  each  of  tee  k  ♦  1  consecutive  problems  has  the  seme  optimal  solution  for  tha 
firm  r  pariods,  timn  this  solution  will  ba  part  of  tha  optimal  solution  to  any  longer 
protease. 


It  Is  also  ponrfbte  to  find  decision  horizons  by  using  another  procedure  when 
the  demand  aaquanoe  k  at  or  near  a  peak  that  k  greater  than  the  available  production 


capacity.  Theorem  16  applies  to  problems  of  length  less  than  or  equal  to  k  and 
Theorem  17  appttes  to  problems  of  length  greater  than  k. 


Theorem  16  OenaUar  the  optimal  aotatioe  to  a  t-period  probiam  t 
whare  t  <  L  If  Xj  -  XMAX|  for  j  -  1.  .  .  .  .  t.  than  period  t  is  a 

t 

fararaar  and  docW—  boric—  for  any  Problem  P* ,  where  t  l  t*  i  T. 

Prooft  If  in  the  optimal  aotart—  to  P‘,  thorn  art  k  consecutive 
periods  with  production  at  the  upper  hound,  than  for  probiam  P* ,  where  t  S 
t*  S  T.  any  i— ada  occurring  aftar  period  t  must  be  aatiafiod  from 
pwdurti—  aftar  period  L  The  moult  b  that  I  «  0  in  any  optimal  solution 
and  the  optimal  production  mpaiea  IX*.  .  .  .  .  X*l  wffl  be  —changed  by 
damanda  occarrtad  aftar  period  t  Tharafora  parted  t  ia  a  dacWon  boric— 
and  ahma  only  t  parioda  of  demand  informati—  am  needed  to  datarmina  it, 
period  t  ia  alee  a  foracaai  horizon,  a 

Theorem  17:  ConaMar  the  optiaaal  sotuti—  to  a  t-pariod  probiam  P* 
wham  t  >  k.  If  Xj  •  XMAXj  for  j  -  t  -  k  ♦  i.  .  .  .  .  u  than  period  t  ia 

a  foraeaet  and  decision  boric—  for  any  Probiam  P* ,  whom  1 1  t*  i  T. 

Pm  aft  If  in  the  optimal  aointi—  to  P\  thorn  am  t  oonaacntivo 
pariodt  haring  prodocti—  at  the  opper  bo— d.  th—  for  problem  P*  wham  t 
d  t*  d  T,  any  damanda  occurring  after  period  t  moat  be  aatiafiod  from 
product!—  aftar  period  t  or  prior  to  period  t  -  k  ♦  L  If  aatiafiod  from 
prodnciti—  prior  to  period  t  -  k  ♦  1,  them  units  will  be  in  inventory  for  at 
knot  k  ♦  1  parioda,  violating  Lemma  1  Therefore  the  damanda  occurring 
aftar  period  t  moat  be  mtisfied  from  product!—  after  period  L  As  in  the 
pmriom  proof,  this  implies  that  I(  *  0  in  any  optimal  soluti—  and  the 
optiaaal  production  acquanco  (X*.  .  .  .  ,  Xjl  win  be  —changed  by  demands 
occurring  after  period  t  Therefore  period  t  is  a  decWoo  boriz—  and  sbtca 
only  t  periods  of  demand  informati—  am  needed  to  determine  it,  period  t  is 


also  a  forecast  horizon.  * 

&  THE  FORWARD  ALGORITHM 

la  this  section,  two  forward  ilioritlmi  arc  described  and  umd  to  solve  Problem 
PT.  Ike  first  dpifta  is  umd  whan  production  capacity  is  constant,  die  second  when 
H  is  not 

We  bathi  by  aakfcw  two  definitions: 

Dafh^Uaa  lit  Let  V(t)  be  the  optimal  value  of  the  t-period 

problem,  t.  wham  t  *  T. 

Daflnitiaw  lb  Let  C(q,t)  be  the  optimal  value  of  the  (t-q*l)-period 
subprobtam,  P(q,t),  that  begins  at  period  q.  ends  at  period  t.  with  die 
rasrrfctions  that  I  *  0.  I  *O,andI>0forj*q.....t-L 

I  j 

The  rcqidicgaent  that  the  P(q,t)  subproMems  have  positive  inventory  in  every 
period  eaoopt  the  laat  is  made  in  order  to  remove  as  much  duplication  of  computational 
effort  as  pomfble  (more  an  this  shortly).  In  effect,  period  t  is  f weed  to  be  the  only 
ragsamation  point  in  Subprobiam  P(q,t).  The  definition  of  Subproblem  P(q.t)  is  similar 
to  that  of  Pierian  and  Klein’s  Capacity  Constrained  Sequences  [23.  Subprobtam  P(q.t) 
can  be  modeled  as  an  integer  linear  program  as  in  Figure  5-L 

Subprobiam  P(q,t)  is  modeled  similarly  to  Problem  PT  as  described  in  Figure 
2-1  except  then  1)  the  only  stockout  variable  is  the  one  in  the  last  period,  S ,  and  2)  a 
penalty  variable  IT  is  included  in  both  the  objective  function  end  in  Constraint  14. 
Let  M  be  a  very  large  number  so  that  if  a  solution  has  Ir  »  0  for  any  r  »  q,  .  .  .  , 
t  -  L  is  equal  to  2*M  and  the  objective  function  attains  t  very  high  value  for  that 


IS 


(10) 


IrH  ♦  Xr  -  Ir  *  dr  for  r«q^..tt-l 

(11) 

I  ♦  X  -  1  ♦  S  -  d 

H  till 

(12) 

Xr  -  XMAXr4r  i  0  f or  r»q^...t  • 

(13) 

M*C|Ir  -  1|  -<Ir  -  1)1  »  0  for  r»q_,t-l 

(14) 

n  *  o 

i 

(15) 

I  -  I  -0 
n  » 

(16) 

Xr,  lr  integer  for  r*q^..,t 

(17) 

St  integer 

(18) 

4  binary  for  r*q^...t 

(19) 

Figure  5*1:  Subproblem  P(q,t) 

Definition  2&  A  solution  to  Subproblem  P(q,t)  with  *  0  for  r  * 
q,  ....  t  -  1  it  aid  to  be  edm/et/b/e.  A  solution  to  subproblem  P(q.t) 
with  tt  least  one  Uf  >  0  for  r  ■  q.  .  .  .  .  t  -  1  it  said  to  be 
Inm/m/uito/B. 
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la  order  to  model  the  absolute  value  condition  in  Constraint  14,  two  new  sets 

of  constraints  are  added: 

Ir  -  (3*r  -  y*r)  *  1  for  r=q^..,t-l  (20) 

y* »  1T  *  0  for  r=q,...,t-l  (21) 

Wo  then  replace  "  j  lf  -  1  j "  in  Constraint  14  by  "(y*  +  fr)n. 

The  nature  of  die  P(q,t)  subproblems  allows  statement  of  several  neccemry 
conditions  given  in  the  previous  section.  In  particular 
L  There  will  not  be  any  stockouts  in  periods  q  through  t  -  1  since  all  the 
inventory  arcs  are  basic  (Lemma  5). 

2.  In  periods  q  ♦  1  through  t.  production  is  either  at  zero  or  at  the  upper 
bound  for  that  period  (Theorem  9). 

1  If  S  >0.  then  X  *  XMAX  (Theorem  10  plus  the  fact  that  X  must  be 
*  q  q  q 

positive  in  order  for  inventory  arc  (q.q+1)  to  be  basic  and  have  positive 
flow). 

4.  If  d  b  XMAX ,  then  there  is  no  solution  to  the  P(q,t)  subproblem,  since 

q  q 

inventory  arc  (q.q+1)  must  be  baric. 

5.  If  the  last  positive  production  occurs  in  period  q'  <  t,  then  d  <  X  S 

j*q’  j  q' 

Z*  d  (Theorem  11  plus  the  fact  that  inventory  arc  (t-l.t)  must  be  basic  with 

n  J 

positive  flow). 


30 

&  In  each  period  ^  else  I,  >  0). 

The  value  V(t)  satisfies  the  following  functional  equation: 

V(t)  *  min  {  V(q-l)  ♦  C(q.t)  )  (22) 

r<q<l 

where  r  *  max  {  t  -  k,  1  }  and  V(0)  3  0. 

We  now  illustrate  how  this  functional  equation  solves  Problem  PT  for  a  small 
nounumerical  example. 

Assume  that  k  *  2  and  T  3  5.  We  begin  by  solving  P1  with  optimal  value 

V(l).  Prom  the  functional  equation,  r  3  1  and 

V(l)  «  V<0>  ♦  cai). 

We  then  proceed  to  solve  the  2-period  problem  P*.  Again,  r  3  1  and 
V(2)  3  min  (  V(l)  ♦  (XU).  V(0)  ♦  C(1.2)  ). 

Notice  that  the  value  of  V(l)  is  already  known  and  that  C(2.2)  is  the  optimal 
solution  to  the  1-period  subproblem  ViU).  The  value  V(l)  +  CX2.2)  immediately  gives 
an  upper  bound  on  the  optimal  solution  to  P2. 

The  3-period  problem  is  solved  next  The  parameter  r  is  still  1  and 
V(3)  3  min  {  V(2)  ♦  CC3.3).  V(l)  ♦  C(2J),  V(0)  ♦  (XU)  h 

The  number  of  possible  completions  to  V(3)  is  3.  which  is  equal  to  k  +  1. 
This  will  always  be  the  maximum  number  of  completions  to  be  evaluated. 

When  V(4)  is  computed,  r  is  now  equal  to  2  and 

V(4)  3  min  {  V(3)  ♦  C(4.4),  V(2)  +  CX3.4),  V(l)  ♦  C<2,4)  ). 
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We  do  not  evaluate  V(0)  ♦  C(1.4)  since  0(1.4)  would  be  the  optimal  value  to 
the  subpreblam  with  poeitivc  eadhif  inventory  in  periods  1.  2.  and  3.  It  it  known  that 
tUa  wffl  never  be  pert  of  an  optimal  sohttkm  to  F4  since  the  maximum  number  of 
consecutive  periods  with  positive  ending  inventory  is'k  which  is  2. 

Finally.  V(5)  is  determined  with  r  *  3  so  that 

V(5)  -  min  {  V(4)  ♦  05.5).  V(3)  ♦  0(4.5).  V(2)  ♦  CttS)  ). 

For  large  values  of  k.  the  computation  of  0(4.0  can  be  quite  formidable  since 
it  is  possible  that  (XMAX^  -  XMIN)  *  2*''  different  production  sequences  would  have 
to  be  enumerated.  However,  in  the  next  section,  a  case  is  treated  in  which  the  P(q,0 
subproblems  can  be  solved  in  polynomial  time. 

One  last  comment  concerns  the  decision  horizon  techniques  and  computational 
savings.  Assume  that  when  we  evaluated  V(4).  we  obtained  e  decision  horizon  by 
Theorem  15  which  covered  periods  1  and  2  (It.,  t*  *  2).  This  means  that  when  we 
evaluate  V(5),  we  need  not  include  periods  1  or  1  Notice  in  the  above  that  when  we 
evaluate  V(5),  we  don't  include  these  periods  anyway,  since  r  ■  3.  The  point-  is  that 
even  if  the  longest  possible  decision  horizon  is  detected  by  Theorem  15.  the  same 
amount  of  effort  must  go  into  evaluating  V(5).  The  value  of  the  horizon  technique 
given  by  Theorem  15.  therefore,  is  not  that  it  saves  computational  effort,  but  that  it 
provides  a  stopping  rule.  That  is.  if  we  only  need  to  determine  the  production  decision 
for  period  1,  we  can  stop  the  forward  algorithm  after  evaluating  V(4)  end  we  need  not 
evaluate  V(5). 

This  is  not  the  case  for  the  decision  horizons  found  by  Theorems  16  and  17. 
If  period  t  *  4  were  found  to  be  s  decision  (and  forecast)  horizon  by  these  methods. 


wu  would  •mam  V(5)  as  a  Mafia  period  problem  consating  only  of  period  5. 

For  tbase  rmsoua,  the  forward  algorithms  in  this  paper  searched  for  decision 
horiaoai  ria  Theoraeas  Id  and  17  first  If  this  was  not  found  to  be  fruitful.  Theorem 
IS  was  mad  to  search  for  decision  horizons. 


5.1.  Computing  (%t)  Whan  Prednctien  Capacity  is  Constant 

When  production  capacity  is  oonstant  from  period  to  period,  then  XMAX.  * 

XMAX  for  all  j.  We  continue  to  assume  that  t  -  q  £  k  or  else  Subproblem  P(q,t) 

wtt  net  ha  part  of  the  optimal  solution.  The  number  of  positive  productions  in 

periods  q  through  t  can  immadktaiy  be  determined  as  follows: 1 

Definition  21:  Let  p  -  I  4  (mod  XMAX). 

rt  J 

DoAnMan  22c  Let  N  «  (Z^d  -  fi/XMAX. 


The  vnhm  of  /  is  the  reeigjndm1  of  the  total  demand  in  periods  q  through  t 
tine  cannot  he  amtafM  by  pmdurtlan  at  XMAX  each  time  period.  Periods  in  which 
pruduetien  is  aa  XMAX  am  called  upper  hounded  production  pmriods.  The  number 
of  appar  hounded  production  periods  in  pvioda  q  through  t  will  be  equal  to  N.  Again 


section  that  if  d^  h  XMAX.  the  Pfq,t)  subproblem  does  not 


Theerum  23c  If  /  *  max  (XMIN,  d  +11.  then  X  *  B,  otherwise  S 

*  0- 

Freed  The  only  positive  production  at  a  quantity  less  than  the 
upper  bound  XMAX  can  occur  in  period  q  (Theorem  9  plus  the  fact  that 


'm  4f  VMM  MlaMaa*  mmh  Mm  Z*  4  d  XMAX  •  (1-4*1 ».  If  mmmpkkm  U  fate*,  ifcaa  mi  X  -  XMAX 

i-4  J  j 

hi  i  •  4,  ....  I.  M  I  «4«M  M  Um  MMilMll  aaflIM  MmM.  mi  cfcack  that  aach  iavaatory  arc  (4,4*1) . 

H*I.U  Mm  pmIDm  Dm. 


each  inventory  arc  is  bask).  It  is  necessary  to  produce  <8  units  in  period  q 
in  order  to  avoid  any  stockouts  in  period  t,  as  long  as  *  XMIN  (Theorem 
11). 

If  fi  <  ♦  1.  however,  inventory  arc  (q.q+1)  will  be  nonbasic  and 

stockout  arc  (S,q)  may  be  baric,  violating  the  definitioo  of  the  P(q,t) 
subproblem.  Therefore  X  *  fi  only  if  fi  >  d  .  If  X  is  not  set  to  0,  the 
first  upper  bounded  production  will  be  in  period  q  and  the  fl  units  will  not 
be  produced  at  all,  but  will  be  stocked  out,  these  stockouts  coming  in  period 
t  e 

The  periods  in  which  to  place  the  N  upper  bounded  productions  must  now  be 
determined.  Their  determination  is  made  possible  by  the  fact  that  the  production  cost 
and  holding  cost  functions  are  constant  over  time. 

The  procedure  for  finding  them  is  to  begin  in  period  q,  determining  X  and  I . 

q  q 

Once  this  is  done,  S  can  immediately  be  determined  (see  above  discussion).  We  then 
move  to  period  t  and  work  beck  to  period  q  ♦  1,  placing  each  XMAX  production  as 
late  as  possible  while  still  keeping  every  inventory  arc  baric.  The  procedure  is 
illustrated  with  a  numerical  example. 

We  wish  to  solve  P(l,4)  where  a  =  5,  p  *  1,  h  =  .3.  s  =  2,  XMAX  =  8.  d|  = 
3,  dj  *  2,  dj  *  6,  and  d4  »  6. 

Period  1  /?  =  17  (mod  8)  =  1  and  N  »  1  Since  fi  <  5  =  XMIN.  X(  =  8.  I(  » 

3,  and  S(  =  0. 

Period  4  1*0  and  S  *  fi  *  1.  Recalculate  d  tobed  -  S  +  I  =  5  (the 

A  A  *  A  A  A  A  ' 


total  inflow  of  production  plus  inventory  into  period  4).  Knee  d4  £ 
XMAX.  set  X4  *  0. 

Mod  3  I#  *  d4  -  X4  *  5  and  S3  =  0.  Now  d,  is  redetermined  to  be  ds  -  S} 
♦  I#  *  1L  Became  d3  >  XMAX,  set  X}  -  XMAX  *  8. 

Piriod  2  Ia  «  d,  -  X#  *  3  and  S2  *  0.  The  demand  d}  will  now  be  replaced 

by  d  -  S  +  I  «  5.  Set  X  =  0  because  d,  S  XMAX. 

2  2  2  2  2 

Finite  Since  I(  *  d2  -  X2.  we  have  found  the  optimal  P(l,4)  subproblem 

solutteo. 

The  solution  to  the  example  follows; 
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3  i 

- 1 
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—  -p  1 
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(  6  1 
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0  i 
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3  1 
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I  0  ; 
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0  I 

0  1 
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l  i  i 
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.  [  — — — 

C(l,4)  «  31.9 

Theorem  24:  The  procedure  described  above  will  find  the  optimal 
solution  to  the  P(q,t)  subproblems  when  production  capacity  is  constant. 

Proof:  If  £  XMAX,  then  Subproblem  P(q,t)  will  have  no 

admissible  solution  since  inventory  arc  (q,q+l)  will  be  nonbasic  with  zero  flow 
in  any  such  subproblem. 


We  first  compute  fi,  the  number  of  units  not  produced  in  upper 
bounded  production  periods,  and  N,  the  number  of  upper  bounded  production 
periods.  The  reasoning  behind  the  placement  of  the  production  (or  stockout) 


X 


Tn  tjrpa  ef  pMm  an  ooaeidand.  ttaa  first  whan  T  S  k  +  1,  the  second 
«mT>  k  ♦  1. 

Mr  a  pnblaa  KiitTlkHw  aohw  for  V(l),  VOX*.  ...  op  to  VCD. 

d turn  the  aMH  a  m  *  Idt*l)/2]  additional  computations  an  performed 

la  am  P*.  Thwef an.  the  Haunt  of  effort  in  soM*  for  VCD  will  be  XT  [t(t+l)]/2 

»•! 

-  OT>*€T>*4T)/l2  or  OCrti 

Pa  a  pnMaa  P1  wH  T  >  k  ♦  1.  we  m>  aom  fa  V(i),  VO),  .  ...  up 
a  vox  but  fOr  P*  wH  t  >  k  ♦  x  WO  only  solve  k  ♦  1  of  the  P(q.t)  subprobktns. 
Thesefae.  we  moat  dtadnpdsh  botwaaa  the  effort  involved  la  solving  P1  through  P**' 
and  Xe  effort  inverted  ia  solving  P**’  fkmagh  PT.  The  effort  involved  in  solving  p1 
Imp  !**'  wB  a  2 £'  [dM)J/2  -  Ok,*12k,*22k+12)/a  The  effort  involved  i» 
aaMat  ^  flnep  P’wib  [<*-DGt-*-2)/2WT-<k*l)].  The  total  effort  is  then 
OClrt  Thh  is  a  vary  Menacing  vault.  U  k  ♦  1  h  T,  the  aoluhoo  algorithm  grows 
«  a  eaMi  tea.  Otherwise,  efar  a  ported  of  cubic  growth,  the  solutioo  algorithm 


U  Caapdai  Cl%t)  Whoa  Pradartiaa  Capacity  is  aet  Cawtaat 

Whaa  da  predacttoa  caparttia  an  aot  coastaat  from  period  to  period,  the 

pverteae  aathad  eaa  aet  be  seed,  lasaad,  the  PCq.t)  subproblems  an  solved  via  a 

breath  sad  hand  tsrtaiqus  Then  an  up  a  (XMAX  -  XMIN)  •  2*^  different 

% 

psahadoa  aequanea  that  way  solve  Subproblem  PCq.0,  each  of  which  must  be 


away  from  consideration. 


Several  procodura  an  employed  a  save  computations  so  that  this  maximum 
production  sequonea  need  aot  be  enumerated.  First,  if  d  *  XMAX .  then 

<1  4 


Of 
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Pitt)  wffi  not  bars  an  admimibte  solution. 

^  for  wUel> x,  *  xma*,  * 

eoaaMerad.  If  a  production  sequence  for  this  reduced  search  has  I.  >  0  for  j  =  q,  .  .  . 
,  t— 1  and  I(  *  0.  then  it  is  an  admissible  solution  to  Subproblem  P(q,t).  If  instead  a 

production  sequence  is  found  for  this  reduced  search  with  I,  >  0  for  j  =  q . t, 

we  check  to  see;  1)  if  every  ending  inventory  I,  for  j  *  q,  .  .  .  ,  t  can  be  reduced 
by  I  units  without  any  of  these  inventories  becoming  zero  (except  l(  itself)  or  negative. 
«d  »  if  X  «.  be  ndwd  by  I(  anils  and  still  btve  greater  Uua  or  equal  to 
XMIN.  If  both  of  these  conditions  are  satisfied,  the  solution  is  modified  by  reducing 
X'  by  I(  units  and  ^  for  j  s  q.  .  .  .  ,t  by  I(  units.  The  result  is  an  admissible 
solution  to  the  P(q,t)  subproblem  with  X  <  XMAX .  The  lowest  cost  admissible 

*  S 

solution  win  be  optimal. 

Theorem  25:  A  reduced  search  on  the  part  of  the  brancb-and-bound 

tree  with  X  *  XMAX  along  with  the  solution  modification  procedure  just 
s  •  ^ 

described  will  generate  the  optimal  solution  to  a  P(q,t)  subproblem. 

Proof:  dearly  die  optimal  solution  to  Subproblem  P(q.t)  will  be 
generated  if  the  solution  has  X  *  XMAX . 

Ajsume  there  exists  an  optimal  solution  V*  to  Subproblem  P(q,t)  with 
first  production  X*  <  XMAX .  It  follows  that  this  solution  win  have  no 

q  q 

stockouts  (Theorem  10).  production  in  periods  q  ♦  1  through  t  at  either  zero 
or  the  upper  bound  (Theorem  9),  and  I*  =  0.  This  solution  will  not  be  in 

•  I 

the  reduced  bnnch-and-bound  tree  because  the  first  production  quantity  is 
not  equal  to  the  upper  bound. 


In  the  reduced  branch-and-bound  tree  however,  (here  win  exist  • 


sofcation  V*  with  XJ  ■  XMAX,  all  other  production  quantities  equal  to  those 
in  V*.  no  stockouts,  and  coding  inventory  levels  that  are  XMAX'  -  X*  higher 
than  those  in  solution  V*. 

If  we  let  A  *  XMAX  -  X*.  solution  V*  can  be  obtained  from 

s  J 

solution  V  by  reducing  the  first  production  quantity  X’  and  every  ending 
inventory  by  A  units.  ■ 

The  definition  of  Subproblem  P(q,t)  can  be  used  to  remove  from  further 

consideration  any  solution  in  which  1^  £  0  for  j  <  L  Solutions  whose  ending 

inventories  are  too  high  can  also  be  removed.  It  is  known  that  in  an  optimal  solution. 

the  ending  inventory  in  a  period  j.  called  1^.  will  be  no  greater  than  the  sum  of  the 

demands  in  periods  j  ♦  1  through  t,  called  IMAX.  or  else  I(  will  be  positive.  However. 

these  ending  inventories  may  eventually  be  reduced  by  I,  units  if  I(  is  positive  in  order 

to  construct  an  admissible  solution.  The  maximum  reduction  that  can  occur  to  these 

inventories  and  still  have  an  admhnibie  solution  is  XMAX  -  XMIN  units,  or  else  X 

*  < 

will  be  tern  than  XMIN.  Furthermore,  for  a  period  j  S  t,  if  I*  =  min  {I  },  then  the 

I  I* 

maximum  reduction  that  can  occur  to  the  inventory  in  period  j  will  be  the 

mmU^UCMAX'-XMINh  If  we  reduce  all  inventories  between  periods  q  and  t  by  I* 
or  more  units,  the  resulting  inventory  level  in  period  t*  will  be  zero  or  negative,  thus 
violating  the  definition  of  Subproblem  P(q,t).  Therefore,  if  I. 

min  {I  a-UCMAX  -XMIN)  >  IMAX,  for  some  j,  the  computation  of  the  solution  to  this 

‘  i  j 

subproblem  can  be  terminated  because  it  is  certain  to  be  bounded. 


The  calculation  of  solutions  to  some  subproblems  can  also  be  eliminated  by 


mteg  east  bounds.  As  V(t)  is  computed  uteg  diffment  velum  of  q  tee  Equation  22). 
an  upper  hound  on  each  new  F(q,t)  subprobkea  is  made  available.  For  example  if  in 
the  aotetion  of  V(5).  it  is  found  that  V(4)  ♦  OS5)  *  2».7.  tten  if  V(3)  *  26.  for 
0(45)  to  be  *  candidate  for  pert  of  the  optimal  solution  to  V<5),  0(45)  meat  be  tam 
than  the  current  bast  nine  of  V(S)  minus  V(3),  or  17.  In  datwminteg  0(4 5).  the 
calculation  of  any  solution  can  be  terminated  that  win  have  a  final  value  greater  than 
or  equal  so  this  upper  bound.  Once  agate,  however,  the  value  of  the  solution  may  be 
overstated  if  it  will  be  reduced  due  to  an  I(  >  0.  The  nettesn  reduction  win  be  by 
(XMAX  -  XMIN)  units.  The  ectual  reduction  win  lower  the  total  production  coat  in 
period  q  and  the  inventory  holding  costs  in  every  period  q  through  t  and  win  be  no 
greater  than  (XMAX  -  XMIN)  •  (p  ♦  h*[t-q*ll).  If  after  making  this  cost 
adjustment,  the  value  of  the  solution,  is  still  greater  than  or  equal  to  the  upper  bound, 
its  calculation  is  terminated  due  to  cost  bounds. 


The  branch-and-bound  technique  is  illustrated  with  an  example.  Assume  we 
wish  to  determine  015)  with  tr  •  5,  p  »  1.  h  ■  .  3,  and  s  *  2.7,  so  that  k  ■  5  and 
XMIN  *  3.  The  demands,  capacities,  end  maximum  inventory  levels  ate  given  in  the 
following  table: 


1  } 

2 

1 

3  : 

4  1 

5  1 

- !• 

■I* 

— —  j 

-  j- 

- 1 

6  I 

9 

1 

4  1 

4  i 

6  i 

9  I 

8 

1 

7  i 

6  I 

7  i 

23  1 

14 

1 

10  1 

6  I 

0  1 

nu 


Besides  die  maximum  inventory  levels,  it  is  also  known  that  7  S  X(  £  9.  X  « 
(O.XMAX }  for  j  -  2.  3.  4.  and  that  X}  «  0  (since  XMAX-  k  dp. 

The  branch-and-bound  tree  for  this  problem  is  constructed  in  Figure  5-1  The 


vataae  on  each  arc  are  the  production  rate*,  ending  inventory,  number  of  stodroun,  and 
coat  id  that  point  ta  the  tret,  respectively.  Several  branches  of  the  tree  are  bounded 

doe  to  iofaadbiHiy  (inventory  going  to  zero  or  stockouts  occurring  before  period  L 

P1M00(L*V*L)  | 

I 

1  j  (9/3/0/14.9) 

/  \ 

/  \ 

2  (8/2/0/28.5)  /  \  (Q/0/6/+00) 

/\  — 

/  \ 

3  (7/5/0/42)  /  \  (0/0/2/+00) 

/\  — 

/  \ 

4  (6/7/0/55.1)  /  \  (0/1/0/42.3) 

I  I 

I  I 

5  (0/1/0/55.4)  "a"  |  |  (0/0/5/55.8)  -b» 


Figaro  5-2:  Branch  end  hound  solution  of  a  P(q,t)  aubprobiem 


Two  ceodldaie  aotadone  are  left  at  the  end  of  the  search.  aohdioo  "a"  with  the 
production  aequenoe  (94.744)  and  aotadoo  "b"  with  producdoo  sequence  (94.7.0.0). 


Sotadon  "a"  bee  aa  extra  enit  hi  inventory  however  and  dun  does  not  aadefy 
the  form  of  Bubproblam  Mq,t).  Thia  aotadoo  h  modified  by  reducing  each  coding 
inveotory  and  the  fkat  period's  producdoo  by  1  unit  with  cost  aavinga  of  15.  A  now 
aotadoo  ia  found  with  prodoctioo  aequenoe  (84.7.6.0).  Thia  aotadoo  baa  a  ooet  of  55.4 
-  15  -  519.  which  ie  opdmaL 


The  optimal  aotadoo,  derived  from  aotadoo  "a"  in  Figure  5-1  ia  given  in  the 


foDewtag  table: 
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Lst  us  no m  look  at  tbo  computational  complexity  question  for  the  variable 
upper  bound  apetahd  lot  doe  model  In  striving  P(q,t),  die  branch-end-bound  tree 
erfll  bare  branches  that  requirint  at  most  that  many  computations.  Again,  the 

value  of  0q,t)  is  not  computed  for  t  -  q  >  k. 


Two  types  of  problama  an  considered,  the  first  when  T  £  k  +  1.  the  second 
when  T  >  k  ♦  L 

For  a  problem  PT  with  T  £  k  +  1,  we  solve  for  V(l).  VO).  .  .  .  ,  up  to  V(T). 
Given  dm  solution  to  P**\  X*  ,  O^'-l)  *  2**'  -  t  -  2  additional  computations  an 
performed  to  solve  P*.  Therefore,  the  amount  of  effort  in  striving  for  V(TX  will  be 
I1,  (21*1  -  t  •  2)  «  2t*j  -  4  -  £Trr+l>J/2  -  2T  or  CX2t). 

For  a  problem  PT  with  T  >  k  ♦  1.  we  spin  strive  for  V(l).  VO).  .  .  .  ,  up 
to  VCD.  but  for  P*  with  t  >  k  ♦  1,  only  k  ♦  1  of  the  P(q,t)  subproblems  an  solved. 
Therefore,  we  mutt  distinguish  between  the  effort  involved  in  striving  P1  through  Pk'f| 
and  the  effort  involved  in  solving  F**’  through  PT.  The  effort  involved  in  solving  P1 
through  F**'  will  be  (Tx  -  t  -  2)  «  2*°  -  4  -  [(k*lXk+2)J/2  -  2<k*I).  The 

i*i 

effort  involved  in  striving  V>*2  through  PT  will  be  [2**2  -  k  -  3MT-<k+l)].  The  total 
effort  is  OCR*).  In  a  fashion  similar  to  the  constant  production  capacity  case,  if  k  ♦  i 
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*  T,  the  solution  algorithm  from  at  an  exponential  rate.  Otherwise,  after  a  period  of 
exponential  growth,  the  solution  algorithm  becomes  asymptotically  linear  in  T. 

t  COMPUTATIONAL  RESULTS  - 

It  has  already  been  shown  that  if  k  +  1  i  T,  then  die  wont  case  solution 

effort  for  the  capeci tiled  lot  sire  model  is  OCT1)  when  production  capacity  is  constant 

and  0(2t)  when  production  capacity  varies  from  period  to  period.  The  average 
performance  of  most  algorithms  is  mually  much  better  than  their  worst  case  bounds  and 
this  is  true  for  the  forward  algorithms  just  danribed.  In  this  section  the  details  of  how 

the  average  computational  performances  were  ^wwIhwI  for  both  die  constant 

production  capacity  and  variable  production  capacity  algorithms  are  given. 

6.1.  Variable  Production  Capacity  Algorithm  Parformanct 

A  group  of  43  test  pro  Warns  wee  created  to  study  the  performance  of  the 
variable  production  capacity  algorithm.  Each  test  problem  wm  50  periods  in  length. 
The  values  of  the  variable  production  cost  per  unit  (p)  end  the  stockout  cost  par  lost 
sale  (s)  wan  kept  fixed  at  $1  and  S3,  respectively.  The  setup  oost  (#)  wm  varied 
between  $40  and  $160.  This  generated  min  liman  positive  production  quantities  (XMINs) 
of  10  units  and  40  units,  respectively.  By  altering  the  holding  cost  per  unit  per  period 
(h),  different  values  of  the  maximum  holding  period  (k)  could  be  generated.  Holding 
costs  of  0.26  (k  *  15).  0.20  (k  *  20),  and  0.16  (k  *  25)  were  used. 

Demands  were  generated  randomly  from  two  uniform  distributions,  U(20.60)  and 
U(0,80).  (We  use  the  convention  that  the  uniform  distribution  U(x.y)  has  minimum 
value  x  end  maximum  value  y.)  Cyclical  demand  patterns  (1  cycle  ■  24  periods)  were 
also  constructed  using  the  sune  ranges  as  the  random  demand  patterns.  Production 
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upper  bounds  ware  p—ad  randomly  from  two  uniform  distributions  as  well,  U(60,80) 
and  U(5830t 

Each  test  problem  was  actually  solved  three  times.  The  first  run  (RUN  I)  used 
both  the  decision  horizon  techniques  and  Subproblem  P(q,t)  fathoming  techniques  (from 
the  previous  section).  The  second  run  (RUN  W  used  only  the  fathoming  techniques. 
The  third  run  (RUN  IH)  used  neither  the  decision  horizon  techniques  nor  the 
fathoming  techniqussL 


The  test  problems  were  solved  on  an  IBM  4341  computer  using  a  forward 
algorithm  written  in  VS  FORTRAN  at  Clarkson  University. 

It  was  found  that  the  marimum  holding  period  (k)  was  the  single  most  factor 
affecting  the  solution  times.  Considering  only  RUN  I,  for  k  *  15  the  average  solution 
time  was  20.6  seconds  (CPU),  for  k  *  20  the  avenge  solution  time  was  2S1.1  seconds, 
and  for  k  ■  25  the  average  solution  time  was  1882.4  seconds. 

Generally,  die  test  problems  with  XMIN  *  40  took  longer  to  solve  than  the 
problems  with  XMIN  *  10.  This  occurred  only  in  RUN  I  and  RUN  II,  which  used 
the  subproblem  fathoming  techniques.  The  differences  in  solution  times  ranged  from 
10%  for  k  »  15  to  33%  for  k  *  25.  This  is  explained  by  the  fact  that  the  higher 
XMIN  value  is  caused  by  the  larger  setup  cost  (<r).  The  larger  setup  cost  results  in 
fewer  setups  and  higher  inventory  levels,  thus  lessening  the  effect  of  fathoming  by 
inventory  lower  bounds.  On  the  other  hand,  this  should  have  increased  the  likelyhood 
of  fathoming  by  inventory  upper  bounds.  The  evidence,  however,  does  not  bear  this 
out 
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As  expected,  the  dedsioo  horizon  technique*  did  not  reduce  the  solution  times 
of  the  test  problem*,  except  in  those  cues  where  horizons  were  found  by  Theorems  16 
and  17.  The  horizon  checking  procedures  sctuaUy  increased  the  solution  times  by  an 
svenfe  of  16* 

The  Subproblem  P(q,t)  fathoming  techniques  had  a  significant  effect  on  the 
solution  times.  Although  RUN  in  wu  not  carried  out  for  k  *  25,  for  k  *  15  and  k 
*  20.  the  average  solution  time  was  81.4%  less  when  using  die  fathoming  techniques 
(RUN  m. 


1  Average  Solution  Tlaes 

1 

in  CPU  Seconds 

k 

XMXN 

!  RUN  I 

RUN  II 

RUN  III 

15 

10 

j  19.489 

19.143 

89.981 

20 

10 

i  220.476 

216.192 

1632.869 

25 

10 

!  1598.864 

1561.156 

— 

15 

40 

i  21.646 

20.788 

90.645 

20 

40 

i  281.667 

275.014 

1629.617 

25 

40 

!  2165.842 

2104.239 

“ 

Table  6-1:  Computational  Results  Variable  Production  Capacities 

Finally,  the  solution  times  for  RUN  II  for  k  *  25  were  regressed  in  order  to 
estimate  the  average  growth  rate.  The  exponential  curve  of  best  fit  was  Y(T)  = 
3.73*(1.41)t  with  r1  *  76.8%,  where  T  is  the  length  of  the  problem  and  Y(T)  is  the 
expected  time  needed  to  solve  a  problem  of  length  T.  measured  in  CPU  milliseconds. 
Although  this  estimate  is  better  than  the  worst  case  bound  of  CX2T),  it  is  still  hopelessly 
exponential 


(I  Cenetent  Production  Capacity  Algorithm  Performance 

In  a  fashion  similar  to  the  above,  a  group  of  24  test  problems  was  created  to 
study  the  performance  of  the  constant  production  capacity  algorithm.  Half  as  many 
last  problems  than  were  used  previously  were  generated  since  the  production  bound  was 
taken  to  be  die  largest  possible  period  demand.  All  other  factors,  including  3  runs  for 

each  test  problem,  were  kept  identical  to  those  in  the  previous  section. 

!  Average  Solution  Tines  in  CPU  Seconds 


k 

XXZM  | 

SUV  I 

SUM  ZZ 

RUM  ZZZ 

15 

10  1 

0.607 

0.595 

0.671 

20 

io  : 

0.789 

0.773 

0.908 

25 

io  : 

0.934 

0.931 

1.098 

15 

40  j 

0.685 

0.694 

0.711 

20 

40  I 

0.884 

0.888 

0.923 

25 

40  ! 

1.037 

1.042 

1.113 

Table  6-1  Computational  Results  Constant  Production  Capacities 

Again,  the  maximum  holding  period  (k)  was  the  single  most  factor  affecting  the 
solution  time.  The  solution  times  were,  however,  much  lower  then  in  the  variable 
capacity  case.  Considering  only  RUN  I,  for  k  *  15  the  average  solution  time  was  0.65 
seconds,  for  k  *  20  the  average  solution  time  was  0.84  seconds,  and  for  k  =  25  the 
average  solution  time  was  0.99  seconds. 

Became  of  the  very  low  solution  times,  any  attempt  to  draw  significant 
conclusions  must  be  prefaced  by  the  remark  that  a  relatively  Urge  pert  of  the  variation 
in  solution  times  could  be  earned  by  the  varying  load  on  the  computer  used.  However, 
the  low  solution  times  did  make  it  possible  to  study  longer  test  problems.  A  new 
group  of  6  test  problems  was  created  that  were  each  100  periods  in  length.  A  random 
demand  pattern  was  used.  Demands  were  generated  from  the  uniform  distribution 


IK20.60).  The  production  capacity  (XMAX)  was  taken  as  60.  The  solution  results  are 
prawn  tad  is  Table  6-3.  Only  die  fathoming  techniques  were  used  in  tbeae  problems 
(no  horizon  checking). 

{  Solution  Tinas 
k  SHIM  i  in  CPU  Seconds 


25 

10 

5 

2.431 

50 

10 

1 

1 

5.839 

100 

10 

1 

1 

9.667 

25 

40 

1 

3.044 

50 

40 

1 

1 

7.305 

100 

40 

1 

11.628 

Table  1-1  Computational  Results;  100  Period  Test  Problems 

Plots  of  the  growth  in  solution  time  ate  given  in  Figures  6-1  for  XMIN  *  10 
sad  in  Figure  6-2  for  XMIN  *  20.  Notice  diet  in  each  case,  die  solution  times  are 
near  equal  for  the  first  k  ♦  1  periods,  after'  which  die  growth  rates  become  linear. 


A  rogwion  was  performed  on  the  two  lest  probtams  with  k  *  100.  It  was 
found  that  tbs  polynomial  curve  of  bast  fit  was  Y(T)  *  0.29T*J  with  rJ  «  97.5%, 
aHghdy  bsttsr  than  the  wont  caw  bound  of  OCT*).  As  before,  T  is  the  length  of  the 
problem  and  Y(T)  is  tbs  expected  time  needed  to  solve  s  problem  of  length  T. 
manmred  in  CPU  milliseconds. 


iMfi  «f  Pnfctan  (Nriofe) 


Figure  6-1:  Solution  Time  Growth:  XMIN  =  10 


